home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / claed0.z / claed0
Text File  |  1996-03-14  |  4KB  |  133 lines

  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAEEEEDDDD0000((((3333FFFF))))                                                          CCCCLLLLAAAAEEEEDDDD0000((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAED0 - the divide and conquer method, CLAED0 computes all eigenvalues
  10.      of a symmetric tridiagonal matrix which is one diagonal block of those
  11.      from reducing a dense or band Hermitian matrix and corresponding
  12.      eigenvectors of the dense or band matrix
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  15.      SUBROUTINE CLAED0( QSIZ, N, D, E, Q, LDQ, QSTORE, LDQS, RWORK, IWORK,
  16.                         INFO )
  17.  
  18.          INTEGER        INFO, LDQ, LDQS, N, QSIZ
  19.  
  20.          INTEGER        IWORK( * )
  21.  
  22.          REAL           D( * ), E( * ), RWORK( * )
  23.  
  24.          COMPLEX        Q( LDQ, * ), QSTORE( LDQS, * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      Using the divide and conquer method, CLAED0 computes all eigenvalues of a
  28.      symmetric tridiagonal matrix which is one diagonal block of those from
  29.      reducing a dense or band Hermitian matrix and corresponding eigenvectors
  30.      of the dense or band matrix.
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      QSIZ   (input) INTEGER
  35.             The dimension of the unitary matrix used to reduce the full matrix
  36.             to tridiagonal form.  QSIZ >= N if ICOMPQ = 1.
  37.  
  38.      N      (input) INTEGER
  39.             The dimension of the symmetric tridiagonal matrix.  N >= 0.
  40.  
  41.      D      (input/output) REAL array, dimension (N)
  42.             On entry, the diagonal elements of the tridiagonal matrix.  On
  43.             exit, the eigenvalues in ascending order.
  44.  
  45.      E      (input/output) REAL array, dimension (N-1)
  46.             On entry, the off-diagonal elements of the tridiagonal matrix.  On
  47.             exit, E has been destroyed.
  48.  
  49.      Q      (input/output) COMPLEX array, dimension (LDQ,N)
  50.             On entry, Q must contain an QSIZ x N matrix whose columns
  51.             unitarily orthonormal. It is a part of the unitary matrix that
  52.             reduces the full dense Hermitian matrix to a (reducible) symmetric
  53.             tridiagonal matrix.
  54.  
  55.      LDQ    (input) INTEGER
  56.             The leading dimension of the array Q.  LDQ >= max(1,N).
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAEEEEDDDD0000((((3333FFFF))))                                                          CCCCLLLLAAAAEEEEDDDD0000((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      IWORK  (workspace) INTEGER array,
  75.             the dimension of IWORK must be at least 6 + 6*N + 5*N*lg N ( lg( N
  76.             ) = smallest integer k such that 2^k >= N )
  77.  
  78.      RWORK  (workspace) REAL array,
  79.             dimension (1 + 3*N + 2*N*lg N + 3*N**2) ( lg( N ) = smallest
  80.             integer k such that 2^k >= N )
  81.  
  82.             QSTORE (workspace) COMPLEX array, dimension (LDQS, N) Used to
  83.             store parts of the eigenvector matrix when the updating matrix
  84.             multiplies take place.
  85.  
  86.      LDQS   (input) INTEGER
  87.             The leading dimension of the array QSTORE.  LDQS >= max(1,N).
  88.  
  89.      INFO   (output) INTEGER
  90.             = 0:  successful exit.
  91.             < 0:  if INFO = -i, the i-th argument had an illegal value.
  92.             > 0:  The algorithm failed to compute an eigenvalue while working
  93.             on the submatrix lying in rows and columns INFO/(N+1) through
  94.             mod(INFO,N+1).
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.